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ABSTRACT 

Context. The interstellar and intra-cluster medium in giant elliptical galaxies and clusters of galaxies is often assumed to be in hy- 
drostatic equilibrium. Numerical simulations, however, show that about 5-30% of the pressure in a cluster is provided by turbulence 
induced by, for example, the central active galactic nucleus and merger activity. 

Aims. We aim to put constraints on the turbulent velocities and the turbulent pressure in the intra-cluster medium of the giant elliptical 
galaxies NGC 5044 and NGC 5813 using XMM-Newton Reflection Grating Spectrometer (RGS) observations. 
Methods. The magnitude of the turbulence is estimated using the Fe XVII lines at 15.01 A, 17.05 A, and 17. 10 A in the RGS spectra. 
At low turbulent velocities, the gas becomes optically thick in the 15.01 A line due to resonant scattering, while the 17 A lines remain 
optically thin. By comparing the (/n.os + /i7.io)/^i5.oi line ratio from RGS with simulated line ratios for different Mach numbers, the 
level of turbulence is constrained. The measurement is, however, limited by the systematic uncertainty in the line ratio for an optically 
thin plasma, which is about 20-30%. 

Results. We find that the (/n.os + /i7.io)//i5.oi line ratio in NGC 5813 is significantly higher than in NGC 5044. This difference can be 
explained by a higher level of turbulence in NGC 5044. The best estimates for the turbulent velocities using resonant scattering and 
upper limits from the line widths, are 320 < V,„,-b < 720 km s"' for NGC 5044 and 140 < K,„rb < 540 km s"' for NGC 5813 at the 
90% confidence limit. 

Conclusions. The high turbulent velocities and the fraction of the turbulent pressure support of >40% in NGC 5044, assuming 
isotropic turbulence, confirm that it is a highly disturbed system, probably due to an off-axis merger. The turbulent pressure support 
in NGC 5813 is more modest at 15-45%. The (/17.05 + /i7.io)/^i5.oi line ratio in an optically thin plasma, calculated using AtomDB 
v2.0.1, is 2cr above the ratio measured in NGC 5044, which cannot be explained by resonant scattering. This shows that the discrep- 
ancies between theoretical, laboratory, and astrophysical data on FeXVII lines need to be reduced to improve the accuracy of the 
determination of turbulent velocities using resonant scattering. 

Key words. X-rays: galaxies - Galaxies: elliptical and lenticular, cD - Galaxies: clusters: intracluster medium - Scattering - Turbulence 
- Atomic data 



E> 1. Introduction by measuring the width of the emission lines. The R eflection 

Grating Spectrometers (RGS, Iden Herder et a n i200l aboard 

The hot X-ray emitting plasma in giant elliptical galaxies and XMM-Newton have the necessary spectral resolution, but be- 

clusters of galaxies is often assumed to be in hydrostatic equi- cause the spectrometers ai-e slitless, the observed lines ai-e 

• librium. In addition to the thermal pressure, there are, how- blurred due to the extended nature of the source. Despite this 

ever, other processes which could contribute significandy to the complication, attempts have been made to put constraints on the 

total pressure of the plasma. One of these sources is turbu- turbulent velocities in groups and clusters of galaxies from the 

lence induced by merger and AGN activity (e.g. [L oeb & Mao width of the lines in RGS spectra jSanders et al.,.2011,) . Upper 

|1994t iRoettiger et ay|1996|; |BoehringereTaD|1993|). Recent es- Umits ranging from 200 km s ' up to 1500 km s"' were found in 

timates from numerical simulations suggest that bubbles blown a sample of 62 clusters, groups and elliptical galaxies, 
by AGN can drive turbulence up to velocities of about 500 km Apart from velocity broadening of lines, there are also indi- 

s (e.g. Briiggen et aU2005HHeinz et al.|2010'). In cluster merg- j-ect methods to constrain turbulent motions in the hot gas, like 

ers, the turbulence can go up to about 1000 km s ' and higher resonant scattering. Although the hot plasma is in general opti- 

(e.g.lDolag et alJ|2005HLau et al. 2009). Typically, turbulence is cafly thin, it can b ecome opticaUy thi ck at the energies of sti'ong 

thought to provide about 5-30% of the pressure support in clus- resonance lines (iGilfanov et alJfl987l) . Since the transition prob- 

ability of such a line is large, it has a high cross section to absorb 

Current X-ray observatories do not have the spectral res- and subsequently re-emit X-ray photons that have the same en- 

olution to measure turbulent velocities in the ICM directly ergy as the line. Along the line of sight toward the core of an 
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elliptical galaxy, resonance-line photons that are emitted in the 
dense central region have a high probability to be scattered out 
of the beam, because the optical depth is large. Due to the con- 
servative nature of resonant scattering, the scattered photons are 
re-distributed to larger projected distances. This effect causes an 
increase of the sur face brightness at the edge of the galaxy (see 
also the review bv lChurazov et al.ll20 10'). The optical depth de- 
pends on the turbulent broadening of resonant lines, especially 
for heavy elements like iron. A comparison of the fluxes of op- 
tically thin and optically thick lines could be used to place con- 
straints on the level of turbulence in galaxies. 

In principle, there are several suitable lines to study resonant 



Table 1. Effective exposure times for the XMM-Newton RGS 
observations after filtering. 



scattering in the X-ray band (see lChurazov et al 



20ia for a list 



of thes e lines emitted by cool and hot clusters). Churazov et alj 
(12004') use the Fe-K complex of lines at 6.7 keV to find a lower 
limit to the turbule nce in the Perseus cluster Inde pendently, us- 
ing the same data, iOastaldello & Molendil (|2004|) reached sim- 
ilar conclusions. The spectral resolution of the current instru- 
ments at these energies is however low, which limits the ac- 
curacy of the measurement. At the moment, the most useful 
lines to study resonant scattering are the 15.01 A, 17.05 A and 
17.10 A lines of FeXVII in the Fe-L complex, whic h are de- 
tectable in relatively cool elliptical galaxies with RGS dXu et al.l 
[2002 ; Werner et al. 2009). The 15.01 A Hne has the largest os- 
cillator strength (/ - 2.31 - 2.73, depending upon the calcu- 
lation) of the lines and is most sensitive to resonant scattering. 
The oscillator strengths of the 17.05 A and 17.10 A are much 
lower, / = 0.12 and / = 10"** respectively. Since all three lines 
are emitted by FeXVII ions, the ratio between the 15 A and 
17 A lines is a good indicator for resonant s catteri ng. In a sam- 
ple of five elliptical galaxies, IWerner et al.l (12009) find that the 
15.01 A line is suppressed in four out of five systems. Since the 
significance of the individual detections is low, 2-4cr per galaxy, 
they could only estimate an upper limit for the turbulent veloc- 
ity in NGC 4636 which is 0.25 times the sound speed at the 
~9Q% confidence level. In these studies, isotropic turbulence is 
assumed. However, these measurements are sensitive to the ra- 
dial direction of the turbulence. If the direction of the motions 
are predominantly radi al or tangential, it would change the re- 
sults (Zhuravleva et al.i i201 lb . 

Recently, long RGS exposures were obtained for the cool 
(kT <1 keV) elliptical galaxies NGC 5044 and NGC 5813. 
These two galaxies have very different morphologies. NGC 5044 
has bubbles which rise buoyantly, but are moved around in all 
directions due to sloshing motions in duced by a recent merger 
(iDavid et al.1120091: iGastaldello et al.l l2009>. NGC 5813, on the 
other hand, has bubbles which rise nearly linearly along the 
NE-SW axis and does not show evidence o f sloshing or re- 
cent merger interactio n dRandall et al.l l20l"ll iT ran et al. 2001; 
lEmsellem etal] l2007h . In the sample of IWerner et al., (,2009 ). 
NGC 5813 shows a hint of a presence of resonant scattering. 
The (/i7.05 + /i7.io)//i5.oi hne ratio is 1.99±0.34 in a 30 ks ob- 
servation with respect to an expected ratio of 1.3, which sug- 
gests that the turbulent velocity is indeed low as is expected 
for a mildly disturbed system. We aim to further constrain the 
(■^17.05 + A7.io)/A5.{)i hne ratio for both galaxies using the deep 
RGS spectra and obtain better limits on the turbulent velocity. 
This would allow us to study the connection between the mor- 
phology of the ICM and the turbulent velocity in two very dif- 
ferent environments. 

In this paper, we use a dist ance of 31.2 Mpc for NGC 5044 
and 32.2 Mpc for NGC 5813 (iTonrv et aLlliOOTI) . At these dis- 
tances, 1' corresponds to 9. 1 kpc for NGC 5044 and 9.4 kpc for 



NGC 5044 



NGC 5813 



ObsID 


Exposure 


ObsID 


Exposure 


0037950101 


20.2 ks 


0302460101 


29.9 ks 


0554680101 


105.3 ks 


0554680201 


54.9 ks 






0554680301 


52.6 ks 


Total 


125.5 ks 


Total 


137.4 ks 



NGC 5813. Metal a bundances are g iven with respect to proto- 
solar abundances bv iLoddeTsI (12003 ) and uncertainties are given 
at the 68% confidence level, unless stated otherwise. 



2. Data analysis 

In our analysis, we use all available RGS data on NGC 5044 and 
NGC 5813 that have exposures larger than 10 ks. Flares were 
discarded by filtering a light curve extracted from the outer re- 
gions of CCD 9. Bins with count rates above 0.25 counts/s were 
discarded. The observations were only modestly affected by soft- 
proton flares. For each observation, less than 20% of the data 
was discarded due to flares. A summary of the observations used 
with their effective exposure time after flare filtering is provided 
in Table [U 

2.1. RGS analysis 

RGS spectra of NGC 5044 and NGC 5813 were extracted us- 
ing the SAS XMM-Newton analysis so ftware version 10.0 (see 
iTamura et al.l 120011 ; Ide Plaa et al.l 12006 , for a detailed descrip- 
tion of the analysis). We use the RGS model backgrounds pro- 
vided by the SAS software as a best estimate for the background 
spectrum and we subtract it from the extracted source spectra. 
Since the suppression of the 15 A line due to resonant scattering 
would be largest in the very centre of the group, we extract the 
spectrum from a 1' wide region (from -0.5' to 0.5' with respect 
to the centre of the group) in the cross-dispersion direction of the 
RGS instrument. 

Since RGS is a slitless spectrometer, we need to take the line 
broadening due to the spatial extent of the source into account. 
The shape of the observed line profile depends on the surface 
brightness distribution of the source in the dispersion direction 
of RGS. In order to model the line broadening, a spatial profile 
is extracted from MOSl data in a region coiTesponding to the 
selected RGS field of view. We put the profile on a wavelength 
grid using the relation AA - 0.138 A AO, where AO is the off- 
set angle in arcminutes with respect to the centre of the source. 
To apply the effect of spatial broadening to the spectral model 
that is calculated during spectral fitting, we use this profile as a 
convolution kernel. The model spectrum is first convolved using 
the spatial profile, before we convolve it with the response ma- 
trix. In the fit, we are able to optimise the width and centroid 
of the observed line profile by changing two parameters of the 
convolution kernel. The multiplication factor s is the width of 
the resulting line profile with respect to the original profile from 
MOS 1 . The line centroid is represented by AA', which is the shift 
of the centroid in Angstrom. The best-fit results of these param- 
eters will be shown in addition to the other spectral parameters 
in the results section of this paper 
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Table 2. Wavelengths and oscillator strengths for all strong 
FeXVII lines in the 15 A to 17 A range. 



Transition A (,Af f f f 



3C 


2s" 2^" 
2s- 2^^ 


'So 


-2s" 


2p^ 


3d 


'Pi 


15.015 


2.31 


2.73 


2.49 


3D 


'So 


-2s2 


2p' 


3d 




15.262 


0.63 


0.61 


0.64 


3F 


2s^ 2p* 
2s- 2p* 


'So 


-2s2 


2p> 


3s 


^Pi 


16.777 


0.11 


0.11 


0.10 


3G 


'So 


-2s2 


2p' 


3s 


'Pi 


17.054 


0.12 


0.13 


0.13 


M2 


2s^ 2p* 


'So 


-2s2 


2p' 


3s 


^P2 


17.097 




5 10-** 


5 lO-** 



Notes. '"^ Wavelengths in A according to NIST; 
lhttp://physics.nist.gov/PhysRefData/ASD/lines_form.ht ml 

Oscillator strengths as provided by Shor^ ( 119791) : NIST 
fhttpV/physics. nist.gov/PhysRefData/ASD/lines_form.html 
I'' Calculated oscillator strengths (D. Liedahl/SPEX) <''> AtomDB 
v2.0.1 dLoch et all2006l) 



2.2. Spectral modelling 

For t he spectral fitting, w e use the SPEX spectral fitting pack- 
age (iKaastra et al.1 IT99d Pl We fit the RGS spectra with a 
single-temperature model of plasma in Collisional lonisation 
Equilibrium (CIE). Since NG C 5044 shows evidence for multi- 
temp erature structure (see e.g. iTamura et aljr2003l : iGrange et alJ 
1201 li) . we also fit its spectrum with a Gaussian D ifferential 
Emission Measure model (gdem, Ide Plaa et al. I l2006h . To mea- 
sure the three lines of Fe XVII at 15 A and 17 A independently of 
the spectral code, these lines have to be removed from the model 
and monochromatic lines have to take their place. In SPEX, this 
can be achieved by ignoring the Fe XVII ion during the calcula- 
tion of the line strengths. In Table |2l we list all strong FeXVII 
lines in the ~15Atol7A range. These lines are modelled with 
delta functions during fitting and also broadened by the spatial 
profile derived from MOSl. Delta lines can be used, because 
the natural width and the Doppler broadening of the lines are 
much smaller than the broadening due to the spatial extent of the 
source. We sum the lines at 17.05 A and 17.10 A, because they 
are not resolved by RGS. 

In the 15-17 A range, there are, of course, also weaker 
Fe XVII lines that contribute to the flux, but which are not mod- 
elled if the Fe XVII ion is ignored. Their contribution turns out to 
be limited. The sixth brightest line in the 15-17 A range, which 
has a wavelength of 15.46 A and is not listed in Table |2] has an 
intensity which is at peak emission only 3.5% of the intensity 
of the line at 15.01 A. Because of the small contribution of these 
weak lines and the contribution of stronger Fe XVII lines outside 
the 15-17 A range, we fix the temperature and Fe abundance to 
the best-fit value that we obtain from a fit where the Fe XVII ion 
is included. 

Fig. □ shows the fluxed RGS spectrum of NGC 5813 in the 
14-18 A band. The model fit does not take into account resonant 
scattering eff'ects. The figure shows that the line at 15 A (3C) 
is overestimated by the model. However, also the other lines 
from Fe XVIII and O VIII are not well fitted. These discrepan- 
cies are most likely due to uncertainties in the atomic data and 
the ionisation balance, which can lead to systematic eiTors up to 
30% in individual line strengths. Clearly, this is an issue for the 
Fe XVII lines as well. Therefore, we discuss the uncertainties in 
the atomic data of Fe XVII extensively in this paper. 



' See |http://vywvy.sron.nl/spexl for more information 



the elliptical galaxies NGC 5044 and NGC 5813 



FeXVII ' ' ' ' FeXVII 

3C 3G+M2 




, I , , , , LJ_j , , , 1 , , L, , , , 1 

14 15 16 17 18 



Wavelength (A) 

Fig. 1. RGS spectrum of NGC 5813 in the 14-18 A band. The 
line properties of the labelled Fe XVII transitions can be found 
in Table |2] 

3. Modelling of resonant scattering 

Modelling the resonant scattering effect that is observed in 
RGS spectra consists of two parts. First, the FeXVII line ratio 
(^17.05 + ^i7.io)/^i5.oi for an optically thin plasma should be de- 
termined. Since we cannot measure this ratio directly from our 
observations due to the yet unknown optical depth of the plasma 
in the 15 A line, it is calculated using atomic data and spectral 
modelling. Secondly, we calculate how this ratio is affected by 
resonant scattering and turbulence in the ICM. 

3.1. The FeXVII line ratios in optically thin plasmas 

The line intensities of the Ne-like FeXVII lines in the 15- 
18 A range (see Table |2]l have been subject to a long stand- 
ing debate in ato mic physics (par t of the debate is summarised 
in, for example, iGu et al.l l2007t) . For many years, the two 
main groups performing laboratory experiments at the National 
Institute of Standards and Technology (NIST) and Laurence 
Livermore National Laboratory (LLNL) were unable to reach 
agr eement on their mea sured Fe XVII line ratios (see, for exam- 
ple lBeiersdorfeiil2003h . On the other side, theorists, using dif- 
ferent methods and codes, also obtained ratios that were neither 
consistent with other theoretical calculations nor with the exper- 
imental data. Only recently, Gillaspv et al. (2011) appear to have 
converged to experimental values that agree with the closest the- 
oretical predictions within 10%. However, given the history of 
this debate, this is likely not the final word on this issue. 

There are several publicly available spectral codes that at- 
tempt to model the line strengths of the Fe XVII lines in astro- 
physical plasmas either based on theoretical calculations, exper- 
imental data, or both. In this work, we compare the results from 
AtomDifl which is integrated in the widely used AP EC model, 
with the SPEX code, which is based on MEKAL (Mew e et aP 
1995h and upda t ed wi th the theoretical Fe XVII calculations by 
Doron & Beh"^(l2002h . The (/i7.05+/i7.io)//i5 01 ratios as a func- 
tion of the temperature for the different codes are plotted in 
Fig.|2l As can be expected from the inconsistencies between the 
experimental atomic data and theoretical predictions, the spec- 
tral codes do not agree. 

- Website: [http://www.atomdb.orgl 
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Fig. 2. Measured and theoretical values for the line ratio 
(^17.05 + ^i7.io)/^i5.oi as a function of temperature. The lines 
indicate the values predicted by the currently available spec- 
tral codes: AtomDB v l.3.1 (grey), SPEX (black) updated with 
lDoron&Behal(l2002h . and AtomDB v2.0.1 (dashed)- The black 
triangle indicates the value found through laboratory measure- 
ments dfieiersdorfer et al. 2004). The grey triangles show mea- 
sur ements of othe r ellipticals performed bv lWerner et al.l (|2009|) 
and lWerner et all ( [2006 ) (M87). For M87, the temperature deter- 
mining the hne ratio is assumed to be the coolest temperature in 
the two-temperature model. The black dots show our measure- 
ments for NGC 5044 and NGC 5813. 



The codes based on older atomic data, AtomDB vl.3.1 
and MEKAL (not shown in Fig. |2]i, both predict a line ratio 
of about 1.0 with a relatively small dependence on tempera- 
ture, dropping from about 1.2 at 0.2 keV to 0.9 at 1.0 keV. 
However, comparisons between stellar flare observations per- 
form ed by the XMM-Newton and Chandra grating spectrometers 
(e.g. iMewe et al. 2001; Ness et al. 2002), laboratory measure- 
men ts (e.g. Beiersd orfer et al. 2004), and theoretical calculations 
(e.g. lDoron & Beha r 2002) suggest that this ratio is probably un- 
derestimated.lDor on & Behan (I2002D found that di-electronic re- 
combination (DR) and resonant excitation (RE) processes are 
dominant in the 0.4-0.6 keV temperature range, while these 
processes were only modelled using scaling laws in MEKAL. 
Because of this finding, SPEX was updated at that time to change 
the predicted Fe XVII line strengths to the new values based on 
these theoretical calculations. The resulting (/17.05 + In.io)/I\5m 
ratio is shown as a black line in Fig.|2] This ratio is considerably 
higher than the results from the MEKAL and AtomDB vl.3.1 
codes. It drops from a ratio of about 1.9 at 0.2 keV to 1.25 at 1.0 
keV. 

The results obtained by iDoron &Behad (l200l are, how- 
ever, not without problems. Although they were right to con- 
sider recombination and ionisation in their model, neither their 
distorted wa ve c alculati ons no r the ionisation balance that they 
adopted (iMazzotta et al.ll998h are state-of-the-a rt for these lines 
or ion s. For example, the single-ion model of iDoron & BehaJ 
(I2OO2I) predicts a value for the (/16.78 + /17.05 + ^i7.io)/^i5.oi ra- 
tio that is always <2 for temperatures below 1 keV. But ac- 
cordi ng to new experimental and theoretical data ( Gilla spv et al.l 
I2OI lb . this ratio should be >2 at these temperatures. The re- 
sults from the 3-ion model (FeXVI, FeXVII, and Fe XVIII) of 



iDoron & BehaH (l2002h are more consistent with the experimen- 
tal values, while they should not be, because the ionisation and 
recombination processes included in this model do not take place 
in the experimental setup. Electron beam ion trap (EBIT) ex- 
periments p redominantly contain io ns from the same kind, i.e. 
Fe XVII. As IDoron & Beh"^ Jlool showed in the 3-ion model, 
the presence of the neighbouring ions (Fe XVI and Fe XVIII), 
as in astrophysical sources, increases the population of the lev- 
els of the 16.78 A, 17.05 A, and 17.10 A lines, thus enhanc- 
ing the (/i6.78 + /17.05 + /i7.io)/^i5.oi ratio with respect to lab- 
oratory measurements. Although this effect is genuine and ex- 
pected in astrophysical plasmas, their predicted values using the 
single-ion model are not consistent with the laboratory values of 
iGillaspv et al.](l20Tll) . 

Recently, AtomDB v2.0.1 was released (iFoster et al .11201 Oh 
containing major updates of the atomic data and new calcula- 
tions of the line formation processes. The Atom DB developers 
chose to use more recent calculations (lBadnelll[2006aib.F I based 
on the code bv lBadnelil (1986). which should be more consistent 
with the experimental data and also includes ionisation and re- 
combination processes, like in Doron & Behar (2002). As shown 
by the dashed line in Fig. |2] the predicted ratio from AtomDB 
v2.0.1 has increased with respect to the other codes. The ratio 
drops from about 2.5 at 0.2 keV to 1.4 at 1.0 keV. 

Comparing the predicted ratios with laboratory experiments 
or astrophysical observations is not straightforward, because it 
is difficult to determine to what extent resonant scattering is act- 
ing in the experimental setup or in the astrophys ical sour c e. The 
tokamak measurements reported by Beiersdor fer et al.l (|2004|) 
are believed not to be affected by resonant scattering. The 
weighted average of their measured data points is shown in 
Fig.|2] The values are consistent within errors with the updated 
spectral codes, but only marginally consistent with MEKAL 
and AtomDB vl.3.1. We also compare these ratios with the 
values measured in a sm all sample of giant elliptical galaxies 
(IWerner et al.ll2009l 12006) . which are shown in grey in Fig. |2] 
If a data point is above the calculated curve, then it is an indi- 
cation of resonant scattering, because the 15 A line can be sup- 
pressed through photons effectively scattering out of the beam 
when pointing to the core of the galaxy. The only data point 
taken next to a galaxy core is "NGC 4636 out". In principle, 
only for this point the 15 A line could be enhanced by photons 
effectively scattered into the beam and lowering the ratio, but 
since the spectrum in this region is dominated by the emission 
near the core of the galaxy where there is no enhancement yet, 
this effect is very small. Otherwise, finding a ratio below the cal- 
culated values is very unlikely. 

The astrophysical data from elliptical galaxies presented 
here appear to favour the ratio found using SPEX above the 
AtomDB v2.0.1 value if we only take the statistical uncertain- 
ties into account. This allows us to favour the numerical value 
obtained by SPEX, but it doe s not m ean that the physical back- 
ground of the lDoron &Behaii(l2002l) model is correct. The older 
MEKAL and AtomDB vl.3.1 underestimate the ratio consider- 
ably, which could be explained by the lower accuracy of the es- 
timations used to model the line formation processes. Although 
the SPEX predictions appear to be more consistent with mea- 
surements, a systematic uncertainty in the line ratio of about 20- 
30% remains when we consider the difference between SPEX 
and AtomDB v2.0.1 (Fig.|2]i and the uncertainty in the oscillator 
strength of the 15.015 A li ne (Table |2ll. This is the same level of 
uncertainty as discussed in lBeiersdorfej (l2003h . 



http ://amdpp .phy s . strath, ac . uk/tamoc/D ATA/ 
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Table 3. Chandra observations of NGC 5044 and NGC 5813. 





NGC 5044 NGC 5813 


ObslD 

Exposure time 


798 9399 
21 ks 84 ks 


5907 9517 
49 ks 100 ks 


Total exposure 


105 ks 149 ks 



3.2. Modelling resonant scattering in the ICM 



Table 4. Optical depths (r) in the FeXVII Unes at 15.01 A and 
17.05 A in the galaxies NGC 5044 and NGC 5813 calculated for 
Mach numbers 0, 0.25, 0.5 and 0.75. 





NGC 5044 


NGC 5813 


Mach 


15.01A 


17.05A 


15.01A 


17.05A 





3.8 


0.2 


6.3 


0.4 


0.25 


1.6 


0.09 


2.6 


0.1 


0.50 


0.8 


0.05 


1.4 


0.08 


0.75 


0.6 


0.03 


0.9 


0.05 



The modelling of resonant scattering in the 15 A line in giant 
ellipticals is performed through Monte Carlo simulations. As an 
input to the simulations, we use deprojected density, temperature 
and Fe abundance profiles obtained with Chandra. The Chandra 
observations of NGC 5044 and NGC 5813 that we use, are 
listed in Table [3] The initial data processing is perforrn ed using 
the standard method described in lVikhlinin et al.l (l2005l) and us- 
ing the latest calibration corrections. Subsequently, we deproject 
the observed spectrum to measure the temperature, density and 
abundance profiles of the plasma. The deprojection technique is 
described in detail by Ch urazo v et al. (2003) and Churazov et aL 
(120081) . In brief, we find the set of 3D spectra in spherical shells 
that fit the observed spectra in concentric annuli best after pro- 
jection. Assuming that the volume emissivity in the outermost 
shell at all energies declines as a power law with radius, we 
model the observed spectra as a linear combination of spectra 
in spherical shells plus the contribution of the outermost layers. 
The matrix that describes the projection of the shells into annuli 
with the contribution of outer layers is inverted and the spectra 
in the shells are calculated by applying the inverted matrix to the 
observed spectra. The resulting s pectra in every sh ell are fitted 
in XSPEC with an APEC model jSmith et al.ll2001h . 

The radial profiles of the deprojected gas parameters in 
NGC 5044 and NGC 5813 are shown in Fig. [3] One can 
see that assuming a fix ed abundance of 0.5 times Solar 
(lAnders & Grevessgll989l) does not change the profiles of the 
density and temperature significantly. The black curves show the 
fits of the density and temperature that we used as an input for 
our Monte Carlo simulations. The number electron density in 
NGC 5044 is described as a yS-model with a normalisation of 
no - 0.06 cm""*, a core radius of rc = 3 kpc and a /? = 0.4. The 
temperature profile (in keV) is parametrised as 



1 + 1.75 



T,{r) = 0.73- 



(1) 



where r is in kpc. The number electron density in NGC 5813 
is also described as a y6-model with a normalisation of no - 
0.15 cm"-', a core radius of rc = 0.5 kpc and a /? = 0.37. The 
temperature is described as 



0.55 




0.75 





0.07 



r > 5 kpc 



-0.09 



(2) 



r < 5 kpc 



3.2.1. Monte Carlo simulations 

In the Monte Carlo simulations, we calculate the optical depth t 
in the centre of the line using 



(Todl, 



(3) 



where I is the distance along the photon propagation direction, 
n, is the number density of ions in the ground state of a given 
transition and cro is the cross section for scattering at the centre 
of a resonant line. This cross-section is given by 



cro 



(4) 



where h is the Planck constant, r^. is the classical electron radius, 
/ the absorption oscillator strength of a given atomic transition 
and AEr, the Doppler width of the line, which is defined as 



AEn = Er 



2kT, 



"^turb 



1/2 



(5) 



Here, £0 is the rest energy of the line, Tg is the electron tempera- 
ture, A is the atomic mass of the ion, nip is the proton mass and c 
is the speed of light. Using the adiabatic sound speed c,, one can 
express t he level of turbulence as a Mach number M = Vturb/cj 
(see, e.g. JChurazov et a02004l) . i.e. 



AEr 



AmpC^ 



(1 + l.4AM^) 



1/2 



(6) 



The calculated optical depths in the considered lines of both 
galaxies are shown in Table |4]for Mach numbers ranging from 
to 0.75. 

Assuming spherically symmetric models of galaxies, we 
then calculate the plasma line emissivities using APEC 
(ISmith et al.l l200lh . Line ener gies and oscillator strengths 
are taken from the AtomDSQ (iLoch et a l. 2006) and NIST 
Database^ The ionisation bala nce (collisional equilibrium) is 
taken from Brvans et al.l (l2009h . Multiple resonant scatterings 
are calculated using a Monte Carlo approach, which treats an 
individual act of resonant scattering in full detail (details of 
these simulations are described in e.g. [C hurazov et al. 2005 
Wer ner et al.ll2009t IZhuravleva et al.ll20iaT201 Ih . In Fig.gl we 
show the resulting ratio of the intensities in the Fe XVII lines, 
(^17.05 + ^i7.io)/A5.oi, for both elliptical galaxies and for Mach 
numbers between and 0.75. One can see the effect of resonant 
scattering: in the core the line intensity is suppressed, while in 
the surrounding regions the intensity of fine increases. 



* http://www.atomdb.org/ 

' http://www.nist.gov/pml/data/asd.cfm 
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Fig. 3. Deprqjected profiles of number electron density, temperature and Fe abundance in NGC5044 (left panel) and NGC5813 
(right panel). Grey points: deprojected profiles with free abundance, black points: deprojected profiles with fixed abundance. Black 
curves: fits of density and temperature used in our calculations. 
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Fig. 4. FeXVII line ratios for NGC 5044 (left) and NGC 5813 (right). The fines give the predicted (/17.05 + /i7.io)//i5.oi line ratio 
for diff'erent Mach numbers (0, 0.25, 0.5, and 0.75) from Monte Carlo simulations a ssuming isotropic tur bulence. Note that these 
profiles have been re-scaled to the optically-thin (/n.os + In.w)/Ii5.oi fine ratio from lDoron & Behan (|2002|) . The data point in each 
plot is the measured line ratio derived from RGS. 



4. Results 

We fit RGS 1 and RGS 2 first order spectra extracted from a 
region 0.5' around the centre of the group with a single temper- 
ature model. In addition, we also fit a multi-temperature gdem 
model to the spectrum of NGC 5044. Unfortunately, the spectral 
resolution deteriorates too quickly to derive meaningful results 
from the outer spatial regions in the cross-dispersion direction. 
The rotation angles of the RGS extra ction regions with re spect 
to the group emission is visualised in lGrange et al.l (1201 ih . The 



best-fit parameters for a single-temperature model and the gdem 
model are listed in Table |5] Due to resonant scattering and the 
systematic uncertainties in the atomic data, the fits are not per- 
fect, but acceptable for our purpose. The multi-temperature fit 
for NGC 5044 is comparable in quality to the single-temperature 
fit and yields parameter values that are consistent with the single- 
temperature results. The width of the lines due to the spatial ex- 
tent of the source are well fitted by the spatial profile extracted 
from EPIC MOS, because the scale factor s (see Section IZTT i is 
consistent with being 1 for both groups, which means that the 
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Table 5. RGS best-fit model results for NGC 5044 and NGC 
5813 for an extraction region of 0.5' around the centre of the 
source. For NGC 5044 we also show the results from a multi- 
temperature gdem model. 



Parameters 


NGC 5044 


NGC 5044 gdem 


NGC 5813 


Y (10™ m-^ 


8.9±0.3 


8.8±0.3 


1.79+0.12 


kT (keV) 


0.784+0.003 


0.795±0.007 


0.647+0.003 


cr 




0.043±0.012 




o 


0.40±0.03 


0.42±0.03 


0.49+0.05 


Ne 


0.81±0.09 


0.72±0.11 


0.32+0.07 


Fe 


0.56±0.02 


0.59±0.03 


0.70+0.05 




1.01±0.03 


1.01+0.03 


0.99+0.03 


A^' (mA) 


-1.1±3.1 


0.2+3.2 


-1.2+2.7 


C-stat / d.o.f. 


959/518 


955/517 


1303/781 



width is the same as the width of the original profile. Also the 
wavelength, which depends on the position of the source, is es- 
timated properly. The wavelength shift AA' is free in the fit and 
the best fit value is consistent with being 0. 

Another reason for the relatively high C-stat values that we 
obtain could be the presence of multi-temperature structure. The 
RGS data of NGC 581 3, however, d o not show significant ev- 
idence for that (Gra nge et al.ll201lh . In NGC 5044, there are 
indications for multi-temper ature structure, but the lower cut- 
off temperature is at 0.4 keV (iGrange et al.ll201 ) h. which means 
that there is no indication for temperatures below 0.4 keV. Fig.|2] 
shows that the expected line ratios flatten out above tempera- 
tures of ~0.5 keV, which just adds an uncertainty of <10% to the 
optically-thin line ratio, which is small compared to the uncer- 
tainties in the atomic data. Since plasma below temperatures of 
~0.5 keV is extremely rare in cool-core clusters and giant ellip- 
ticals (ISanders & F abian 2011) and we do not see evidence for 
OVII emission lines indicating temperatures below ~0.4 keV, 
we will use the results from the single-temperature model. 

To determine the (/n.os + /i7.i())//i5.()i line ratio, we use the 
best-fit model shown in Table |5] and ignore FeXVII in the spec- 
tral calculation. This removes the strongest FeXVII lines, which 
are listed in Table|2] from the model. We therefore replace these 
five lines with four delta functions. Since we cannot resolve the 
17.05 A and 17.10 A lines, they are modelled by a single delta 
function at the average wavelength of the two lines. Then, we 
fit the spectrum again with the temperature and Fe abundance 
fixed to their best-fit values, because these two parameters are 
likely to introduce bias due to the lack of Fe XVII lines in the 
spectral model (See Sect. 12.21 1. Other parameters, like the nor- 
malisation, the oxygen and neon abundances, the delta lines, and 
the spatial broadening parameters are allowed to optimise. From 
these fits, we obtain a (In 05 +/17 io)/^i5 01 line ratio of 1 .25±0.08 
for NGC 5044 and 1.87±0.08 for NGC 5813. Also the ratio 
from the multi-temperature fit of NGC 5044 yields 1.25±0.08, 
which shows that the line ratio does not depend on the multi- 
temperature structure. The ratios measured in these two galaxies 
are significantly different from each other, which suggests that 
the resonant scattering effect in NGC 58 1 3 is much stronger than 
in NGC 5044. 

The comparison of the observed ratios with the line ratios 
from an optically-thin plasma calculated by the different spec- 
tral codes yields a very interesting result. In Fig.|2] the observed 
ratio for NGC 5044 appears to be below the predicted ratio from 
AtomDB v2.0.1 at a ~ 2cr confidence level considering only 
statistical errors. Since the (/17.05 + In.\o)/I\5.m ratio can only 



Table 6. Estimates for the turbulent velocities in NGC 5044 and 
NGC 5813. 



NGC 5044 NGC 5813 





Fred." 


M 


Vtuib" 


Fred.* 


0.00 


1.76 


0.00 





2.94 


0.25 115 


1.53 


0.25 


104 


2.13 


0.50 230 


1.42 


0.50 


208 


1.75 


0.75 345 


1.37 


0.75 


311 


1.61 


00 00 


1.27 


CO 


CO 


1.30 


Observed 


1.25+0.08 


Observed 


1.87+0.08 


Lower Limit 


> 320kms-' 


Lower Limit 


> 140kms-' 


Upper Limit 


< 770kms-' 


Upper Limit ' 


< 540 km S-' 



Notes. Equivalent turbulent velocity estimate in km s ' using the 
Mach number and the adiabatic sound speed, which is 460 km s"' for 
NGC 5044 and 415 km s"' for NGC 5813. Emission weighted aver- 
age of the predicted (/17.05 + /i7.io)/^i5.oi ratio over the RGS extraction 
region. 90% upper limits on the turbulence derived from the line 
width in RGS. 



increase due to resonant scattering in the core of an elliptical 
galaxy, we cannot explain the low value for NGC 5044 if the 
AtomDB v2.0.1 ratio is assumed. The v alue of 1 .25 is, however , 
consistent with the predicted ratio using lDoron & Beh aij (l2002|). 
The el liptical ga laxy NGC 4649 in the sample of |Werner et^ 
(2009) and M87 (i Werner et a l.'2006') also show a ratio of ~ 1 .25- 
1.30, but with a larger uncertainty. It should be noted that one 
of the points, NGC 4636 out, was derived from a RGS spec- 
trum just outside the galaxy centre where the optical depth of 
the 15 A line is expected to be smaller than in the centre. 

Since the (/17.05 + /i7.io)/^i5.oi ratio for an optically- 
thin plasma is uncerta in from the atomic data side and the 
iDoron & Behad (l2002h results are problematic (Section 13.1b . 
we are left without a good reference ratio to base our turbu- 
lence estimates on. The best we can do is to assume that the 
hot gas in NGC 5044 is subject to such high levels of turbu- 
lence that the gas is effectively optically thin and estimate the 
turbulence in NGC 5813 relative to that ratio. Since the shape 
of th e (/17.05 + /i7.io) //is.oi ratio as a function of temperature 
from iDoron & Behaij (|2002) is very similar to the shape of the 
AtomDB v2.0.1 curve (see Fig.|2|, and because it is also con- 
sistent with the ratio from NGC 5044, we assume for this paper 
that Doron & Behar ( 2002) indeed describe the Fe XVII ratios in 
an optically-thin plasma correctly and that the gas in NGC 5044 
is optic ally thin. At the same tim e, we are aware of the problems 
with the lDoron & BehaJ (l2002h result. Therefore, we discuss the 
systematic uncertainties extensively. 

4. 1 . Turbulent velocity estimates 

The (/17.05 + /i7.io)/A5.oi line ratios from the fits can be directly 
compared to the predicted ratios from our Monte Carlo simula- 
tions. The observed ratios and the Monte Carlo results are plotted 
in Fig.m Note that the simulation results have been re-scaled to 
be relative to the optically-thin line ratio predicted using SPEX. 
The horizontal error bar on our data point indicates the RGS ex- 
traction region. From the simulations, it is evident that the pre- 
dicted ratio varies substantially within this region. In order to 
make a fair comparison, we choose to calculate the emission- 
weighted average of the predicted line ratio within the -0.5'- 
0.5' field of view. We use the count rates observed in the cross- 
dispersion direction of RGS in this interval as weighting factors 
for the calculation. The resulting ratios are listed in Table|6] 
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Using the emission weighted ratios from the simulation and 
the observed ratios in Table |6] we can constrain the level of the 
turbulent velocity in the inner 0.5' (r < 5 kpc) of the galaxy. For 
NGC 5044, this is more difficult than for NGC 5813, because 
of our assumption that the gas in NGC 5044 is optically thin. 
At these low optical depths, the resonant scattering method loses 
sensitivity in the velocity estimate. Therefore we can only derive 
a lower limit for the turbulence. If we only consider statistical 
errors, the 90% confidence lower limit is 320 km s ' for NGC 
5044. For NGC 5813, the 15 Aline is significantly suppressed by 
resonant scattering with respect to the SPEX prediction, which 
allows a better estimate of the turbulence. To obtain the veloci- 
ties in km s"', we interpolate linearly between the simulated val- 
ues in Table|6l Since the errors due to the linear interpolation are 
smaller than the systematic errors in the atomic data, we do not 
gain accuracy by performing more Monte Carlo simulations for 
other Mach numbers. Our best estimate for Vturb in NGC 5813, 
without taking into account any systematic error, is 175+22 km 
s"', which means a statistical precision of about 13%. 

From the discussion in Sect. [3j however, we know that we 
also need to take into account the systematic uncertainty in the 
optically-thin line ratios. For the calculation above, we used the 
line ratio predicted by SPEX, which appears to be a reason- 
able estimate. But, the typical difference between SPEX and 
AtomDB v2.0.1 is about 15%, which also adds to the uncer- 
tainty. Using the AtomDB v2.0.1 ratio would not give a phys- 
ical result for NGC 5044, but we can calculate it for NGC 5813. 
Assuming an optically-thin ratio of 1.513, the turbulent veloc- 
ity would be 310 km s"' . The difference between the velocities 
based on the AtomDB v2.0.1 and the SPEX ratio is a factor of 
two, which means that deriving the turbulent velocity in this way 
is very sensitive to the assumed optically-thin line ratio. 

Although the optically-thin line ratio is sensitive to the 
plasma temperature, we do not expect that multi-temperature 
structure in the gas adds to the systematic uncertainty. Fig. |2] 
shows that the line ratio has its largest decline between 0.2 and 
0.5 keV. From 0.5 to 1.0 keV, the ratio flattens. The tempera- 
tures that we measure for NGC 5044 and NGC 5813 are in the 
range 0.6-0.8 keV, where the ratio is fairly constant. The ratio 
could be biased if there were temperatures present below 0.5 
keV. However, we do not detect significant O VII emission in 
the spectrum, which would indicate the presence of those low 
temperatures. It is therefore safe to assume that the Fe XVII line 
ratios are not biased due to a significant contribution from low 
temperatures and that the systematic uncertainty in the derived 
velocities due to multi-temperature structure is negligible. 

Another method of estimating the turbulence from high- 
resolution spectra is by measuring the line width. Although the 
line broadening due to the spatial extent of the source makes an 
accurate determination of the turbulent velocity impossible, it is 
possible to obtain an upper limit of the turbulent velocity with 
RGS. In the fit, we both allow the turbulent velocity in the CIE 
model and the scale factor s that controls the width of the spatial 
distribution to be free. Since the spatial distribution is different 
from being Gaussian, it does not correlate significantly with the 
Gaussian velocity broadening, which allows us to obtain upper 
limits with a good accuracy. During the error calculation, also 
the scale factor s is allowed to optimise, which means that un- 
certainty is included in the upper limit. The results are shown on 
the last row of Table|6] For NGC 5044, the 90% upper limit con- 
strains the turbulent velocity to the range between 320 < Vturb < 
770 km s-i. The upper Hmit for NGC 5813 of < 540 km s"' is 
consistent with the estimated value from resonant scattering. 



An important question to understand the structure of the ICM 
in groups and clusters of galaxies is what fraction of the pressure 
in the ICM is due to turbulence. This fraction is equivalent to the 
ratio of the turbulent energy (eturb) to the thermal energy (etherm), 
which can be easily calculate d from the Mach number of the 
turbulence dWerner et al.ll2009h : 

= (7) 

^therm 2 

where y is the adiabatic index for an ideal gas (y - 5/3). From 
the lower limit derived for NGC 5044, we derive that the tur- 
bulent energy can range from 40% of the thermal energy up to 
a factor of 2 times the thermal energy. The turbulent pressure 
clearly is a component which should be considered for NGC 
5044. The turbulent energy fraction for NGC 5813 is signif- 
icantly lower, between 15% and 45% if we take into account 
the systematic uncertainty in the initial line ratio. Since we only 
probe the inner 0.5' of both galaxies with RGS and the optical 
depth in the 15 A line decreases rapidly above r > 10 kpc, the 
values derived above are valid for the inner regions within r < 
5-10 kpc. 

4.2. Chandra images of NGC 5044 and NGC 5813 

Since we find a significant difference between the measured tur- 
bulent velocities in NGC 5044 and NGC 5813, it is interesting 
to look for the origins of the turbulence by looking at the spatial 
structure of the groups. They both show X-ray cavities blown 
by the central AGN in their core, which is a possible source of 
turbulence. In order to visualise these cavities, we used the 4 
Chandra ACIS observations listed in Table [3] to create an im- 
a ge where th e average ICM emission is subtracted (see also e.g. 
IDavid et al.ll20C^ [Randall et alJl20riL for NGC 5044 and NGC 
5813 respectively). In all observations, the source was centred 
on ACIS-S3 of the ACIS-S detector. For both sources, the two 
individual event files from the standard pipeline products were 
merged to create a single event file containing the maximum pos- 
sible exposure. Then, we extracted radial surface brightness pro- 
files within a 2.5' radius from the central AGN in the 0.5 - 7 
keV band. The maximum extraction radius of 2.5' falls within 
the CCD 7 boundaries in all observations. The derived surface 
brightness profiles were fitted using an empirical model consist- 
ing of two King profiles and a constant. The full expression for 
this model /(r) is: 



/('•) = 











ft '■'] 
















i >-cl \ 









with normalisations A^i and N2, core radii and r^a, and expo- 
nents Pi and 132 for both King profiles. A constant C was added 
to account for the background emission. The expected number 
of counts from the best-fit double-King model was subsequently 
subtracted from the observed number of counts in the original 
image for every pixel. The residual images that we obtained were 
smoothed with a Gaussian using a width (cr) of 0.98". The result 
for both galaxies is shown in Fig.|5] 

5. Discussion 

Using deep XMM-Newton RGS data of NGC 5044 and NGC 
5813, we have found a significant difference between the (/17.05 + 
^i7.io)/^i5.oi hne ratios in these galaxy cores. Our result shows 



8 



J. de Plaa et al.: Estimating turbulent velocities in the elliptical galaxies NGC 5044 and NGC 5813 




that the FeXVII line at 15 A in NGC 5813 is significantly af- 
fected by resonant scattering, confirming the 2cr hint found by 
IWerner et alj ( 12009 ) that is based on a shorter 30 ks RGS expo- 
sure. The (/i7.05 + /i7.io)/^i5.oi line ratio in NGC 5044, on the 
other hand, is found to be consistent with the line ratio predicted 
by the CI E model in SPEX, wh ich is in agreement with the con- 
clusion of lTamura et alJ (l2003h . based on a 20 ks RGS exposure, 
that resonance scattering effects in NGC 5044 are small. The 
difference between the line ratios in these galaxies is best ex- 
plained by resonant scattering in the 15 A line of NGC 5813, 
which is strong due to the relatively modest turbulent velocities 
(~0.4M) in the core of this object. Other effects that could af- 
fect the measured line ratio, like multi-temperature structure, do 
not appear to play a role. Interestingly, we do not find strong 
evidence for multi-temp erature structu r e in NGC 50 44, contrary 
to the results found by 'Tamur a et alJ (I2003D and Gran ge et alJ 
(l20Tll) . This difference originates from our relatively narrow 
cross-dispersion selection that we need to detect resonant scat- 
tering. We extract the RGS spectru ni from a V w i de reg ion in 
the cross-dispe rsion direction, whil e iTamura et alj (l2003h use a 
width of 2' and'Grang e et al.l (1201 1') a width of 5'. The origin of 
the multi-temperature signal reported in these previous works is 
therefore most likely the temperature gradient in the inner few 
arcminutes of the source. 



5.1. Systematics in the optically-thin line ratio 

The systematic uncertainty in the unscattered optically-thin 
(^17.05 + ^i7.io)/A5.oi ratio is the most important limitation for 
measuring turbulent velocities in clusters and elliptical galax- 
ies. For NGC 5813, we find a difference of a factor of two in 
the derived turbulent velocity when the optically-thin line ratio 
is increased by only 15%. This clearly adds a large systematic 
uncertainty to the derived turbulence. The result can only be im- 
proved when the uncertainty in the optically-thin line ratio de- 
creases, which is difficult due to the problems with the atomic 
data of FeXVII. The estimated uncertainty on individual line 
strengths produced by spectral codes is about 20-30%. 



RGS observations of elliptical galaxies like NGC 5044 also 
help to constrain the optically-thin line ratio. Although we cur- 
rently lack the statistics to definitively rule out certain model 
calculations, a larger sample of deep observations of these ob- 
jects will help to constrain the upper limits for this ratio. As 
shown in Section 13.1! laboratory measurements and theoreti- 
cal calculations could eventually improve our knowledge about 
the optically-thin line ratio as soon as they yield consistent re- 
sults. There are signs that theore tical and laboratory results are 
converging (lGillasDvetal.ll2011h . but this still needs confirma- 
tion. Considering the remaining uncertainties in individual line 
strengths, it is clearly necessary to invest in theoretical and ex- 
perimental research to improve the atomic data b efore new mis- 
sions like Astro-H are launched dTakahashi et al.. .2010.) . This 
work shows that accurate atomic data is important to perform 
robust measurements of astrophysical sources. 



5.2. Line profiles and velocity broadening 

Upper limits on the turbulent ve locity were es t imated for NGC 
5044 and NGC 5813 before bv ISanders et"al] dMT ) based on 
the oldest data sets with unfiltered exposures of 43 and 73 ks, 
respectively. They find upper limits of 1500 km s ' for NGC 
5044 and 1400 km s ' for NGC 5813 at the 90% confidence 
level. Our upper limits are significantly smaller. For NGC 5044, 
we find an upper limit based on line broadening of <770 km 
s"', which is a factor of ~1.95 improvement with respect to 
Sanders et alJ (l201 ll) . This is consistent with the expected im- 
provement (~ VS) due to the increased statistics. For NGC 5813, 
the improvement is much larger than expected based on statis- 
tics. The upper limit drops from <1400 km s"' to <540 km s"' 
(factor ~2.5), which is much better than the improvement of the 
statistics (factor V2). We note that the method of correc ting for 
the spatial broadening is diff'erent in'Sande rs et al.l (1201 ih . there- 
fore the results may not be fully comparable. Although our upper 
limit for NGC 5813 is lower than expected, it is still consistent 
with the independent measurement using resonant scattering. 



9 



J. de Plaa et al.: Estimating turbulent velocities in the elliptical galaxies NGC 5044 and NGC 5813 



We must also realise that the velocity broadening measure- 
ments of elliptical galaxies with RGS are sensitive to systematic 
effects. The line profile of a line in RGS does not only depend 
on the velocity broadening, but also on the line strength as a 
function of radius in the elliptical galaxy. Since the line strength 
of a certain line depends on the temperature, density, element 
abundance, and resonant scattering properties, the line profile 
observed with RGS can be different for each line. In order to in- 
clude the bias due to the different line widths in the upper limit, 
we allow the width of the broadening profile derived from MOS 
(s) to optimise. The derived upper limits on the turbulence thus 
give a rough estimate of what level of turbulence would still be 
consistent within the average of the fine widths of the lines in the 
spectrum. 

In addition, resonant scattering would also change the pro- 
file of the 15 A line. The core of the line would be suppressed 
and because resonant scattering is a conservative process, the 
line wings should be enhanced. Unfortunately, line blends pre- 
vent us from studying the wings of the 15 A line in more detail. 
Since the strongest resonant scattering effect is the suppression 
of the core of the line, our method of using delta lines convolved 
with the spatial line profile give us a good estimate of the line 
ratios. The other ICM properties that determine the line profile, 
like temperature, density, and abundance do not influence our 
estimate of the (/n.os +/i7.io)/-^i5.oi ratio, because these lines are 
all emitted by Fe XVII and they should have a similar profile in 
the observed temperature range, apart from resonant scattering 
effects. 



5.3. Origin of turbulence 

The RGS spectra of NGC 5044 and NGC 5813 were extracted 
from the bright inner regions (r < 5 kpc) of both galaxies. The 
optical depth in the 15 A line declines rapidly beyond a radius 
of ~10 kpc. It is therefore reasonable to assume that most of 
the resonant scattering signal originates from the inner <10 kpc 
of the galaxy. Strictly speaking, any differential gas motion can 
contribute to the line broadening and affect the optical depth. 
Therefore, our estimate of the turbulence and the turbulent en- 
ergy implicitly includes all types of different motions, provided 
that they happen on spatial scales where the accumulated opti- 
cal depth is >1. If one assumes normal turbulence (e.g. with a 
Kolmogorov spectrum) most of the energy is found at the largest 
scales. Only eddies with a characteristic size comparable to or 
smaller than the size of the optically thick co re affect the am- 
phtud e of the resonant scattering significantly (IZhuravleva et al.l 
I20T1I) . 

Both NGC 5044 and NGC 5813 show AGN activity in the 
form of bubbles, therefore the considerable difference in the 
measured turbulence may be difficult to explain by AGN induced 
turbulence alone. The AGN activity in NGC 5813 has pr oduced 
bubbles and relatively strong shocks (iRandall et al.ll201 1). but 
apparently only modest differential gas motions and turbulence 
along the line of sight compared to NGC 5044. The inner 10 kpc 
of NGC 5044, on the o ther hand also shows e vidence of recent 
~ 10^ yr AGN activity (iGastaldeflo et al.ll2009l) . but no evidence 
of strong shocks like in NGC 5813. Therefore, the AGN activity 
in NGC 5813 appears to be stronger than in NGC 5044, which 
is opposite to the conclusion that one would reach based on the 
limits on the turbulence that we obtained. AGN activity alone 
cannot explain the difference in turbulence that we detect be- 
tween NGC 5044 and NGC 5813. 



Another source of turbulence would be cold fronts induced 
by an off axis merger Chandra, XMM-Newton, and optical 
data show evidenc e that NGC 5044 h a d a recent me rger with 
a smaU sub group (iMendel et aLll2008t iGastaldeUo et al.li2009l: 
iDavid et al.l2009l) . The Chandra images of NGC 5044 and NGC 
5813 in Fig.|5]indeed show a striking difference. While the bub- 
bles in NGC 5813 appear to rise in a straight line along the NE- 
SW axis, the bubbles in NGC 5044 rise in many directions. Most 
likely, the recent merger in NGC 5044 has induced a sloshing, 
or swirling, motion in its ICM, causing the bubbles to rise in 
a rotating medium jDavid et al.l 12009). Along the edges of the 
cold fronts associated with the merger, instabilities cause turbu- 
lent motions that may be detected using resonant scattering. A 
potential problem with thi s interpretation is the fact t hat nu mer- 
ical sim ulation s , Uke e.g. lAscasibar & Markevitc B (12006 ) and 
Roediger et al.l (1201 ih show that for a sloshing cluster the ve- 
locity field in the centre of the hot gas has a much smaller am- 
plitude than the velocities in the outskirts. Since we are probing 
mainly the inner ~10 kpc of the galaxy, we should not see the 
high velocities associated with the merger On the other hand, 
these simulations were made for larger systems and different 
merger geometries. We cannot exclude the possibility that in 
NGC 5044 the gas in the centre is disturbed by the sloshing 
motions. NG C 5813 does not ap pear to be disrupted by a re- 
cent merger (Randall e t al.ll20l"l !). Most likely differential gas 
motions, including turbulence, due to the gas sloshing are re- 
sponsible for the observed small optical depth of the FeXVII 
line at 15 A. 

The derived values for the turbulent energy fraction in NGC 
5813 (15^5%), assuming isotropic turbulence, ar e compara- 
ble to estimates from n umerical simulations (e.g. Dola g et aP 
120051: iNagai et al.ll2007h . which range between 5-30% of the 
thermal energy. The turbulent velocities around 100-200 km s ' 
could very well be induced by the bubbles blown by the AGN. 
Heinz et al. (2010) show that bubbles can generate turbulence as 
high as ~500 km s"' in clusters of galaxies. Due to the reduced 
scale of an elliptical galaxy with respect to a cluster, an AGN 
generated velocity of 100-200 km s ' seems reasonable. NGC 
5044, on the other hand, shows a very high turbulent energy frac- 
tion >40%. The ICM of NGC 5044 appears to be very disturbed. 
Due to the low mass of NGC 5044, a merger is the most likely 
explanation for the high level of turbulence. 



5.4. Future prospects 

Estimating turbulent velocities through resonant scattering and 
line broadening, like presented in this paper, is at the limit 
of the current instrument capabilities. Although fitting high- 
resolution spectra from new micro-calorimeter instruments, like 
SXS aboard Aifro-// (Takahashi et al. 2010), will be a challenge 
for the current atomic databases and spectral codes, the excel- 
lent spectral resolution of SXS will enable us to measure tur- 
bulent velocities directly from the line widths. Since the spec- 
tral band of SXS ranges from 0.3-12 keV, it will also allow the 
study of resonant scattering in other strong resonance lines, for 
example in the Fe-K region around 6 keV (see the review by 
IChurazov et al ]l2010h . The combined effort of improving atomic 
data and flying new instruments like SXS aboard Astro-H will 
enable a more comprehensive study of turbulence in elliptical 
galaxies and clusters of galaxies. The high spectral resolution 
will not only allow measurements of the amplitude of the mo- 
tions, but also their anisotropy and spatial scales. 
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6. Conclusions 

We analyse deep XMM-Newton RGS spectra of NGC 5044 
and NGC 5813. Following a detailed analysis of the (/n.os + 
^i7.io)/^i5.oi line ratio, we find that: 

- The (/i7.05 + /i7.io)//i5.oi ratio in NGC 5813 (1.87±0.08) is 
significantly higher than in NGC 5044 (1.25±0.08). We in- 
terpret this difference to be line suppression of the Fe XVII 
line at 15 A due to resonant scattering, which strongly de- 
pends on the turbulent velocities in the ICM. 

- The core of NGC 5044 shows a (/i7.o5-h/i7.io)//i5 .oi line ratio 
that is consistent with the theoretical ratio by Do ron & Behad 
(|2002 >). It is, however, 2cr below the AtomDB v2.0.1 ratio, 
which cannot be explained by res onant scattering . Although 
the distorted wave calculations bv lDoron & Beha r (2002) are 
not state-of-the-art, it is consistent with our ratio in NGC 
5044. Therefore, we assume this optically-thin line ratio to 
estimate the turbulent velocities in NGC 5044 and NGC 
5813. 

- Using the measured (/17.05 + /i7.io)/^i5.oi ratio, we estimate 
the level of turbulence within r < 5-10 kpc of the centre 
of these galaxies assuming isotropic turbulence. The core 
of NGC 5813 shows moderate turbulent velocities between 
140< Vtuib <540 km s"' (90% confidence) with a best esti- 
mate of ~ 175 km s"'. The ratio of the turbulent energy over 
the thermal energy, which is the same as the fraction of the 
turbulent pressure support, yields a turbulent energy fraction 
between ~15^5%, which is comparable to estimates of the 
turbulent pressure support from numerical simulations. The 
lower and upper limit for the turbulence in NGC 5044 yields 
a velocity range of 320< Vtui-b <770 km s"\ which means a 
turbulent energy fraction of > 40%. 
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